Mechanisms of pattern formation during T cell adhesion 
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T cells form intriguing patterns during adhesion to antigen-presenting cells. The patterns at the 
cell-cell contact zone are composed of two types of domains, which either contain short TCR/MHCp 
receptor-ligand complexes or the longer LFA-l/ICAM-1 complexes. The final pattern consists of 
a central TCR/MHCp domain surrounded by a ring-shaped LFA-l/ICAM-1 domain, while the 
characteristic pattern formed at intermediate times is inverted with TCR/MHCp complexes at the 
periphery of the contact zone and LFA-l/ICAM-1 complexes in the center. In this article, we 
present a statistical-mechanical model of cell adhesion and propose a novel mechanism for the T 
cell pattern formation. Our mechanism for the formation of the intermediate inverted pattern is 
based (i) on the initial nucleation of numerous TCR/MHCp microdomains, and (ii) on the diffu- 
sion of free receptors and ligands into the contact zone. Due to this inward diffusion, TCR/MHCp 
microdomains at the rim of the contact zone grow faster and form an intermediate peripheral ring 
for sufficiently large TCR/MHCp concentrations. In agreement with experiments, we find that the 
formation of the final pattern with a central TCR/MHCp domain requires active cytoskeletal trans- 
port processes. Without active transport, the intermediate inverted pattern seems to be metastable 
in our model, which might explain patterns observed during natural killer (NK) cell adhesion. At 
smaller TCR/MHCp complex concentrations, we observe a different regime of pattern formation 
with intermediate multifocal TCR/MHCp patterns which resemble experimental patterns found 
during thymozyte adhesion. 

PACS numbers: 



I. INTRODUCTION 

T lymphozytes mediate immune responses by adher- 
ing to cells which display foreign peptide fragments on 
their surface. These peptide fragments are presented 
by MHC molecules on the cell surfaces, and recognized 
by the highly specific T cell receptors (TCR). Cytotoxic 
T cells directly kill virally infected cells after adhesion, 
while helper T cells help to activate other lymphozytes, 
e.g. B cells which produce and secrete antibodies. 

At the contact zone of a T cell with an 
antigen-presenting cell (APC), bound receptor- 
ligand pairs are arranged in characteristic 
supramolecu lar patterns, th e 'immunological 
synapse' iMonks et al. 19981 iGrakoui et al. 19991 
iKrummel et al. 20fKTTpcn^eret^ L_2QlH^ for reviews see 
Ivan der Merwe et al. 200(1 iPnstin and Cooner 200f 
Delon and Germain 2000. iBromlev et al. 2 001. 

Dustin et al. 200 ll IWiilfing et al. 2002j . A current 
focus in T cell biology is to understand the formation 
of these patterns and their role for T cell activation 
|K.-H. Lee et al. 2002j . The final pattern of an ad- 
hering T cell consists of a central domain with bound 
TCR/MHCp complexes, surrounded by a ring-shaped 
domain in which the integrin receptors LFA-1 of the 
T cell are bound to their ligands ICAM-1 of the APC. 
Intriguingly, the characteristic intermediate pattern 
formed earlier during T cell adhesion is inverted, with a 
TCR/MHCp ring surrounding a central LFA-1/ICAM- 
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1 domain in the contact zone [Grakoui et al. 199*91 
l.lohnson et al. 200(1 IK.-H. Lee et a,l. 2002j. This nat- 
tern inversion has been first observed for T cells adhering 
to a supported lipid bilayer with MHCp and ICAM-1, 
|Crakoui et al. 19991 l.lohnson et_gL_2(2QCj [ . more rece ntly 
also in a cell conjugate system |K.-H. Lee et al. 2002j . In 
the case of natural kille r (NK) cells, the inverted synapse s 
seem to be stable [Davis et al. 19991 iFassett et al. 200l) . 
More recently, thymozytes (immature T cells), have been 
found to form multifocal synapses with several nearly 
circular clust ers of TCR/MHC-peptide complexes in the 
contact zone [Hailman et al. 20021 iRichie et al. 2002j . 

As has been proposed by several groups, the lateral 
segregation or phase separation in the immunological 
synapse probably is caused by the length difference 
between receptor /ligan d co mplexes iDavis et al. 1996 , 
Shaw and Dustin 1997, Ivan der Merwe et al. 2000, 
Qi et al. 200ll IS.-.T. Lee et al. 20021 IWeikl et al. 2002a . 
Burroughs and Wiilfing 20021 IChen 2003 , 

CoornbsetHdlf ! Bound TCR/MHCp complexes in- 
duce a membrane separation of about 15 nm, while 
LFA-l /ICAM-1 complexes have a larger extension of 
40 nm [Dustin and Cooper~2 000] . In general, the lateral 
phase behavior of the membranes is also affected by 
the concentrations of the c omplexes | Weikl et al. ' 2002a. 
IWeikl and Lipowskv 200~ll IWeikl et al. 2002b| . Lateral 
phase separation only occurs if the complex concentra- 
tions exceed a critical threshold. An additional driving 
force for phase separation comes f rom large glycopro - 
teins such as CD43 and CD45 [Weikl et al. 2002rj| . 
These glycoproteins ha ve a length of 40 nm and more 
[Shaw and Dustin 1 997] and thus form a steric barrier 
for TCR/MHCp binding. 

To understand the time-dependent pattern evolution, 



it is necessary to consider the dynamics of the phase sep- 
aration process in the cell geometry. The first theoreti- 
cal model based on a Landau-Ginzburg free energy cou- 
pled to a set of reaction-diffusion equatio n has been pro- 
posed by Qi, Grove s, and Chakraborty |Qi et al. 200ll 
IS. -J. Lee et al. 2002j . Qi et al. obtain circularly sym- 
metric adhesion patterns exhibiting the characteristic do- 
main inversion and conclude that the entire T cell pattern 
evolution may be caused by self-assembly. The formation 
of the intermediate inverted synapse is suggested to be 
formed in a pivot mechanism in which the formation of 
LFA-l/ICAM-1 complexes in the center of the cell con- 
tact zone causes close apposition in the periphery allow- 
ing TCR/MHCp binding. 

Here, we consider a statistical-mechanical model for 
the cell adhesion dynamics in which the T cell and APC 
membranes are discretized into small patches of linear 
size a ~ 70 nm. In this model, the configurational energy 
of the membranes depends on the numbers of receptors, 
ligands, and glycoproteins present in each patch, and the 
separation of apposing patches in the cell contact zone. 
The adhesion dynamics is studied with Monte Carlo sim- 
ulations, which include the whole range of fluctuations in 
membrane shape and composition. Based on the pat- 
terns observed in the simulation, we propose a novel 
mechanism for the formation of the intermediate inverted 
synapse. We observe an initial formation of many small 
TCR/MHCp microclusters throughout the contact zone. 
Subsequently, TCR/MHPp clusters at the periphery of 
the contact zone grow faster due to the diffusion of free 
receptors and ligands into this zone. For sufficiently large 
TCR/MHCp concentrations, these peripheral clusters co- 
alesce into a closed ring, surrounding a central domain of 
LFA-l/ICAM-1 complexes. At smaller TCR/MHCp con- 
centrations, we observe a different dynamic regime with 
characteristic multifocal intermediate patterns consisting 
of several circular TCR/MHPp domains, which resemble 
patterns observed during thymozyte adhesion. 

Our proposed mechanism for the formation of the in- 
termediate inverted T cell synapse is a self-assembly 
mechanism based on TCR/MHCp microcluster nucle- 
ation and the diffusion of free receptors and ligands into 
the contact zone. However, we find that the final, 'ma- 
ture' T cell synapse only arises in the presence of active 
cytoskelctal transport of TCRs to the center of the con- 
tact zone. This seems to be in agreement with experimen- 
tal findings. Cytochalasin D, an inhibitor of actin-based 
transp ort, inhibits also th e central TCR/MHCp move- 
ment [Grakoui et al. 199 9], An active actin/myosin- 
based transpor t of receptors into the c ontact zone has 
been observed |Wulfing and Dav is 1998], while the gly- 
coprotein CD43 is actively moved out of th e contact zone 
|Allenspach et al. 200ll iDelon et al. 20011] . Without ac- 
tive transport, the intermediate TCR/MHCp ring seems 
to be metastable and persists for an hour and more, ac- 
cording to our simulations. This might explain the in- 
verted synapse of natural killer (NK) cells which con- 
sists of a peripheral ring of short HLA-C/KIR complexes 
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FIG. 1: Cartoon of a T cell membrane (top) adhering to 
an APC membrane (bottom). The T cell membrane con- 
tains the T cell receptor TCR (green) and the receptor LFA-1 
(red). The APC membrane contains the corresponding lig- 
ands MHCp (green) and ICAM-1 (red). Both membranes 
contain repulsive glycoproteins (grey). Because of the differ- 
ent lengths of bound TCR/MHCp complexes, LFA-l/ICAM-1 
complexes, and glycoproteins, the membrane phase separates 
into domains. 



surrounding a central domain containing the longer LFA- 
l/ICAM-1 complexes. The formation of the inverted NK 
synapse is not inhibited by ATP depletion or disruption 
of the cytosk eleton and thus appears to be cause d by 
self-assembly |Davis et al. 1999tlFassett et al. 2001") . 



II. MODEL 

In this section, we describe our theoretical model for 
the interaction of a T cell with an antigen-presenting cell 
(APC). We consider two apposing membranes. The first 
membrane represents the T cell and contains the recep- 
tors TCR and LFA-1. The second membrane represents 
the APC and contains the ligands MHCp and ICAM-1. 
We use the terms 'receptors' and 'ligands' here with re- 
spect to the T cell: Adhesion molecules anchored in the 
T cell membrane are called 'receptors', and those in the 
APC membrane are 'ligands'. Protruding glycoproteins 
are embedded in both membranes, forming a steric bar- 
rier for the formation of the short TCR/MHCp complexes 
(see Fig. Pi . 

To mimick the adhesion geometry of the cells, we di- 
vide the membranes into a contact zone and a surround- 
ing region in which the membranes do not interact. The 
receptors can diffuse in the whole T cell membrane, but 
interact with the ligands of the APC membrane only 
within the contact zone of the two membranes. For sim- 
plicity, we avoid the problem of modeling the full cell 
shape, and assume here that the contact zone has an es- 
sentially circular shape and a constant area on the time 
scales considered here, see Fig. This contact zone is 
thought to be established in fast initial adhesion events 
after first cell contact. Experimental pictures of adhering 
T cells show tha t the contact zone fu lly develops in less 
than 30 seconds [Grakoui et al. 1999j . 

To characterize the membrane conformations, we par- 
tition both membranes into quadratic patches with linear 
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extension a. 1 The local composition of the T cell mem- 
brane is then described by the numbers nf of TCRs, nf 
of LFA-1, and nf* of glycoproteins in each membrane 
patch i. Correspondingly, the composition of the APC 
membrane is given by the numbers nf 1 of MHCp, nf of 
ICAM-1, and nf a of glycoproteins in all patches. Within 
the contact zone, the local separation between two appos- 
ing membrane patches of the two cells is denoted by U. 

The elastic energy of the membranes in the contact 
zone is dominated by the bending energy and by lateral 
tension and can be written as 



^ciW = [(K/2a 2 )(A d k) 2 + (a/2)(V d k) 



(1) 



Here, k — kik^/^i + ^2) denotes the effective bending 
rigidity of the two membranes with rigidities n\ and K2, 
and a is a lateral tension. For simplicity, the effective 
bending rigidity is taken to be independent of the local 
membrane composition. The term 
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is the total curvature of the membrane separation field U 
at site i, and 
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describes the local area increase of the curved mem- 
branes with respect to the reference x-y plane given by 
h — l x ,y — 0. The elastic energy Q dominates the fluc- 
tuations of the membrane separation in the contact zone, 
whereas the overall cell shape is also affected by the elas- 
ticity of the cytoskeleton which is coupled to the mem- 
branes. In the simulatio ns, wc use the dimensionless sep- 
aration field z = {l/a)^J 'k/ '(fc^T), and choose the value 
z = 1 to correspond to a length of 20 nm, which results 
in the relation a — 20^ k/^bT) nm for the linear patch 
size. For the typical rigidities n\ — K2 — 25ksT of the 
two biomembranes, the effective rigidity k has the value 
12.5/cbT, and the linear patch size is a ~ 70 nm. For 
the lateral tension, we choose the value a = 0.1k/ a 2 ~ 
2 • IQ-^N/m. 

The overall configurational energy of the membranes 
in the contact zone is the sum of the elastic energy Q 
and the interaction energies of receptors, ligands, and 
glycoproteins: 

H{1, n} = H cl {l} + [^hw + min(nf , nf)Vru{k) 

i 

+ wm(nf ) ni)V u (l i ) + (nf* + nf") V G (k)] (2) 

Here, Vtm(^) and V\a{U) are the attractive interac- 
tion potentials of TCR/MHCp and LFA-1 /ICAM-1 com- 
plexes, respectively, while Vg{U) is the repulsive in- 
teraction potential of the glycoproteins. The term 



1 More precisely, we are discretizing the reference plane shown in 
Fig.EJinto a square lattice with lattice constant a, which results 
in partitioning both membranes into quadratic patches. 




FIG. 2: 'Cell' adhesion geometry: The circular contact zone 
is surrounded by a nonadhering membrane ring. Receptors, 
ligands, and glycoprotein diffuse around in the whole mem- 
brane, but interact with the apposing membrane only within 
the contact zone. 



min(nf, nf 1 ) denotes the minimum of the numbers of 
TCR and MHCp molecules at site i. This minimum 
is equivalent to the number of interacting TCR/MHCp 
pairs in the apposing patches at site i. 

The receptor-complexes can only form if the membrane 
separation is in an appropriate range. The length of 
the TCR/MHCp complexes is about 15 nm, while the 
LFA-l/ICAM-1 complexes have a length of about 40 nm. 
Since the membrane within a patch is 'rough' due to the 
thermal fluctuations on length scales smaller than the lin- 
ear extension a ~ 70 nm of the patches, we assume that 
the complexes can form if the separation of two appos- 
ing patches does not deviate more than 5 nm from the 
lengths ztm and zli- The interaction potential of TCR 
and MHCp is then given by 

Vtm — Utm for 10 nm < It < 20 nm, 

= otherwise (3) 

and the interaction potential of ICAM-1 and LFA-1 is 

Vlj = Uli for 35 nm < U < 45 nm, 

= otherwise (4) 

where Utm < is the binding energy of a TCR/MHCp 
complex, and Uli < the binding energy of LFA- 
1/ICAM-l. As noted above, the potential width of 
10 nm effectively takes into account small-scale fluctu- 
ations within patches. Thus, this width does not re- 
sult from the atomic interaction potentials of receptor 
and ligand molecules, which should have a significantly 
smaller range. Similarly, the binding energies of the 
receptor/ligand complexes should be seen as effective 
binding energies which can be used to adjust the 2- 
dimensional equilibrium constants of the TCR/MHCp 
and LFA-l/ICAM-1 complexes. The 2-dimensional equi- 
librium constants are approximately given by Ktm — 

a 2 e -U TM /(k B T) and „ a 2 e -U LI /(k B T) (gee Ap _ 

pendix) . 

The repulsive glycoproteins protruding from both 
membranes vary in size. However, many of these pro- 
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FIG. 3: Time sequence of Monte Carlo conformations of the contact zone for the effective binding energies Utm = — 6.5/cbT 
of TCR/MHCp complexes and Uli = -3k B T of LFA-l/ICAM-1 complexes. The overall concentrations of TCR, ICAM-1, 
LFA-1, and glycoproteins in each of the membranes is OA/a 2 ~ 80 molecules//im 2 for a linear patch size a ~ 70 nm, and the 
concentration of MHCp is 0.1/a 2 ~ 20 molecules//im 2 . Membrane patches with bound TCR/MHCp complexes are shown in 
green, patches with bound LFA-l/ICAM-1 complexes in red. 



teins have a length comparable to the length of the LFA- 
l/ICAM-1 complexes. These glycoproteins do not in- 
hibit the binding of ICAM-1 and LFA-1, but impose 
a steric barrier for the formation of TCR/MHCp com- 
plexes. They are characterized here by the potential 

V a = U G (l - l G f for I < l G , 

= otherwise (5) 

with Uq = 10k/ci 2 and Iq — 40 nm. This potential 
results from the fact that a membrane patch of size a 
containing a glycoprotein has to bend around this protein 
to achieve an overall patch separation smaller than the 
length of the glycoprotein. 

In the following, the radius of the circular contact zone 
is chosen to be 45a, and the nonadhering membrane sur- 
rounding the contact zone is a ring of width 55a. As 
boundary condition at the rim of the contact zone, the 
membrane separation is fixed at a value of 100 nm, which 
is significantly larger than the length of the TCR-MHCp 
and LFA-l/ICAM-1 complexes and the glycoproteins. 

III. ADHESION DYNAMICS IN THE ABSENCE 
OF CYTOSKELETAL TRANSPORT PROCESSES 

We first consider the pattern formation in the absence 
of active forces which transport molecules in or out of 
the contact zone. The lateral motion of receptors, lig- 
ands, and glycoproteins within the membranes is then 
purely diffusive. In our discretized membranes, the dif- 
fusive motion of the macromolccules is modeled as a 
hopping process between neighboring membrane patches. 
Each receptor, ligand, or glycoprotein in a certain mem- 
brane patch can hop to one of the four nearest neighbor 



patches during a single time step. The hopping processes 
of macromolecules located in the nonadhering membrane 
region do not change the cell interaction energy, eq. 
However, within the contact zone, the attempted hop- 
ping of a macromolecule may change the free energy. 
According to the standard Metropolis criterion, the hop- 
ping attempt is always accepted if it does not increase 
the free energy, but is only accepted with probability 
exp(— AF , /(fcsT)) if it leads to a free energy increase 
AF. 2 During a time step, we also attempt to shift the 
separation U between apposing membrane patches in the 
contact zone by d ■ £[—1,1] where d is the step width 
10 nm, and ^[ — 1, 1] is a random number between —1 and 
1. 

A single Monte Carlo step roughly corresponds to 1 ms 
of real time. This time estimate can be derived from the 
2-dimensional diffusion law (x 2 ) = ADt and the typical 
diffusion constant D ~ 1 /im 2 /s for membrane-anchored 
macromolecules. In a single Monte Carlo step, a free 
receptor, free ligand, or a glycoprotein moves a distance 
a to a neighboring membrane patch, which corresponds 
to a diffusion time t = a 2 /(AD) ~ 1 ms for a = 70 nm. 
On the length scale of our patches, the diffusive motion of 
the macromolecules is slower than the relaxation o f the 
membrane separation [Brochard and Lennon 1975| and 
hence defines the time scale. 

These dynamic rules and the free energy given in 



2 We reject moves in which bound ligands or receptors hop from 
one binding partner to another in a single time step. For these 
moves, the free energy difference would be zero. Thus, the ac- 
tual free energy barrier for the unbinding process of the lig- 
and / receptor complex would not be captured. 
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FIG. 4: Time sequence of Monte Carlo conformations of the contact zone for the effective binding energies Utm = — 5.5/cbT 
of TCR/MHCp complexes, Uli = — 4fc_gT of LFA-l/ICAM-1 complexes, and the same molecular concentrations as in Fig. [3] 
Membrane patches with bound TCR/MHCp complexes are shown in green, patches with bound LFA-l/ICAM-1 complexes in 
red. 



eq. J5J specify a stochastic adhesion process. Here, 
we study the adhesion process with Monte Carlo 
simulations. Taking averages over many indepen- 
dent Monte Carlo runs gives a nu merical solution of 
the corresponding master e quation [van Kampen 19921 
iBinder and Heermann 19 92], The stochastic process 
captures the fluctuations in the membrane separation 
and describes the diffusive motion of the receptors, lig- 
ands, and glycoprotein on a single-molecule level, which 
is essential for the mechanisms of pattern formation con- 
sidered in this article. As initial conformation, we choose 
the separation profile I = l a + cr 4 where r is the distance 
from the center of the contact zone, l is 45 nm, and 
c > is chosen so that the separation at the rim of the 
contact zone with radius r — 45a is 100 nm (boundary 
condition). This initial separation in the contact zone is 
larger than 45 nm, and thus beyond the interaction range 
of receptors, ligands, and glycoproteins. Initially, these 
molecules are taken to be randomly distributed within 
the whole membrane. 

We systematically study the adhesion dynamics for 
various concentrations of the receptors, ligands, and gly- 
coproteins and for various effective binding energies, or 
2d equilibrium constants, of the TCR/MHCp and LFA- 
l/ICAM-1 complexes. Since the length difference of the 
complexes leads to phase separation at the molecular 
concentrations considered here, the two types of recep- 
tor/ligand complexes have to 'compete' for the contact 
zone. In general, the overall area of TCR/MHCp do- 
mains in the contact zone increases with the concentra- 
tions of TCR and MHCp molecules and with the effective 
binding energy Utm- However, if the molecular concen- 
trations or the binding energy are too small, TCR/MHCp 
domains do not form, and the contact zone contains only 
bound LFA-l/ICAM-1 complexes. At molecular concen- 
trations and binding energies where TCR/MHCp and 



LFA-l/ICAM-1 domains coexist, we observe two differ- 
ent regimes for the dynamics with clearly distinct pat- 
terns of TCR/MHCp domains at intermediate times. 
The pattern evolution roughly depends on the overall 
area of TCR/MHCp domains after initial relaxation. 

Regime 1: If the overall area of TCR/MHCp domains 
is relatively large, we observe a characteristic ring-shaped 
TCR/MHCp domain at intermediate times, surrounding 
a central domain of LFA-l/ICAM-1 complexes. A typ- 
ical example for the pattern evolution in this regime is 
presented in Fig. [21 The first Monte Carlo snapshots 
of the contact zone show the formation of many small 
TCR/MHCp microclusters. At later times, the clusters 
close to the rim of the contact zone grow faster, and form 
an intermediate peripheral TCR/MHCp ring. The faster 
growth of the clusters close to the rim is caused by the 
diffusion of unbound TCR and MHCp molecules from 
the nonadhering membrane into the contact zone. Fi- 
nally, the ring breaks to form a single large TCR/MHCp 
domain. 

Regime 2: For smaller TCR or MHCp concentrations, 
or smaller effective binding energy, we observe charac- 
teristic multifocal TCR/MHCp patterns at intermediate 
times. A typical example is shown in Fig.0] Initially, we 
observe again the nucleation of many TCR/MHCp micro- 
clusters throughout the contact zone. However, the over- 
all area of TCR/MHCp domains now is not large enough 
for the formation of a TCR/MHCp ring. Instead, mi- 
croclusters in the whole contact zone grow and coalesce, 
which leads to multifocal intermediates and finally again 
to a single TCR/MHCp domain. 

To distinguish the two dynamic regimes systematically, 
we consider a peripheral ring of the contact zone with 
distances r > 35a from the center, and divide this ring 
into 100 equal segments. For each Monte Carlo pattern 
obtained during adhesion, we determine the fraction Y 
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of ring segments which contain bound TCR/MHCp com- 
plexes. A fully closed peripheral TCR/MHCp ring corre- 
sponds to a ring occupation Y = 100 %. We find that the 
crossover between the two dynamic regimes can be ap- 
propriately described by a maximum ring occupation of 
ymax — go % attained during adhesion. A pattern evolu- 
tion with y max < 80 % typically has multifocal interme- 
diates as in Fig. 0] (Regime 2), while pattern evolutions 
with y max > 80 exhibit the inverted synapse of T cells 
with peripheral TCR/MHCp ring, see Fig. [21 (Regime 1). 

The diagram at the top of Fig. 03 shows how the dy- 
namic regimes for pattern formation depend on the effec- 
tive binding energies Utm and Uli of the TCR/MHCp 
and LFA-l/ICAM-1 complexes. These binding energies 
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FIG. 5: Dynamic regimes for T cell adhesion. The con- 
centrations of TCR, LFA-1, and ICAM-1 are OA/a 2 ~ 80 
molecules/pm 2 and the concentration of MHCp is 0.1/a 2 ~ 20 
molecules/^im 2 . In the top diagram, the glycoprotein concen- 
tration in each of the membranes is Xg = OA/a 2 . In the bot- 
tom diagram, the binding energy Uli of LFA-l/ICAM-1 com- 
plexes has the value 3fcsT. The black diamonds in the figure 
represent data points obtained from Monte Carlo simulations. 
At large values of the binding energy Utm of the TCR/MHCp 
complexes, we observe a peripheral TCR/MHCp ring at inter- 
mediated times as in Fig. |H] (Regime 1). At medium values of 
Utm, multifocal patterns as in Fig.^Jare obtained at interme- 
diate times (Regime 2). At small values of Utm, TCR/MHCp 
domains in the contact zone do not form. The threshold for 
the formation of TCR/MHCp domains and the crossover be- 
tween the two dynamic regimes depend on the binding energy 
Uli of LFA-l/ICAM-1 complexes and the glycoprotein con- 
centration Xg in both membranes. 





FIG. 6: Two possible patterns with a single TCR/MHCp do- 
main, shown in green. In the left pattern, the TCR/MHCp 
domain only has boundaries with the red LFA-l/ICAM-1 do- 
main. In the right pattern, the TCR/MHCp domain is in 
contact with the rim of the contact zone. 



are proportional to the logarithm of the 'ideal' 2d equi- 
librium constants of the complexes, see Appendix. An 
increase in Utm in general leads to more TCR/MHCp 
complexes in the contact zone, while an increase in Uli 
leads to the binding of more LFA-l/ICAM-1 complexes. 
We observe three different scenarios: (i) At small values 
of Utm, TCR/MHCp domains do not form at all in the 
contact zone which then is completely occupied by LFA- 
l/ICAM-1 complexes. TCR/MHCp domains only form 
above a threshold value for Utm- This threshold value 
increases with Uli- (u) At large values of Utm, we ob- 
serve Regime 1 of pattern formation with the characteris- 
tic peripheral ring of TCR/MHCp complexes as in Fig. [3] 
(iii) At intermediate values Utm, we find the patterns of 
Regime 2 with characteristic multifocal intermediates as 
in Fig.0] The crossover value of Utm separating Regime 
1 and Regime 2 increases with Uli- The intermediate 
TCR/MHCp ring of Regime 1 only forms if sufficiently 
large numbers of TCR/MHCp complexes are present in 
the contact zone. Instead of varying the effective bind- 
ing energies Utm and Uli, the numbers of bound recep- 
tor/ligand complexes in the contact zone could also be 
changed by varying the overall concentrations of the re- 
ceptors and ligands, with similar effects on the pattern 
formation. 

The diagram at the bottom of Fig. \5\ shows the effect 
of the glycoprotein concentration Xq on the adhesion 
dynamics. The length of the glycoproteins is compatible 
with the length of the LFA-l/ICAM-1 complexes. Hence, 
the glycoproteins can enter the red LFA-l/ICAM-1 do- 
mains in the contact zone, but are excluded from the 
green TCR/MHCp domains. The accessible membrane 
area for the glycoproteins increases with the fraction of 
LFA-l/ICAM-1 domains in the contact zone, and so does 
the entropy of the glycoprotein distribution. Therefore, 
an increase in the overall glycoprotein concentrations 
leads to a larger fraction of red LFA-l/ICAM-1 domains 
in the contact zone, and thus has a similar effect as in- 
creasing the binding energy Uli of the LFA-l/ICAM-1 
complexes. 

In both dynamic regimes of pattern formation, the co- 
alescence of clusters finally leads to a single TCR/MHCp 
domain in our model. In the absence of active transport 
processes, we always observe that the final TCR/MHCp 
domain is in contact with the rim of the contact zone, see 
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FIG. 7: Pattern formation with active transport of TCRs towards the center of the contact zone. Membrane patches with 
bound TCR/MHCp complexes are shown in green, patches with LFA-l/ICAM-1 complexes in red. Molecular concentrations 
and binding energies are the same as in Fig. [3] (ring regime). (Top) At zero force, the intermediate TCR/MHCp pattern is 
stable for 30 minutes and more. In the final equilibrium pattern, both types of domains are in contact with the rim of the 
adhesion region, see section III. (Middle) At the force F — O.OlfcsT/a, the final equilibrium state is the target-shaped mature 
synapse of T cells. This state is already established within 30 minutes. (Bottom) At the 10-fold stronger force F = 0.1 ksT/a, 
the final target-shaped pattern already forms within 5 to 10 minutes. An intermediate pattern with a TCR/MHCp ring appears 
around 30 seconds after initial contact. 



Figs. |2| and This behavior can be understood from the 
line tensions at the domain boundaries and at the rim of 
the contact zone. The line tension between TCR/MHCp 
and LFA-l/ICAM-1 domains, A, is an energy per unit 
length of the domain boundaries, and is mainly caused by 
the elastic energy of the membanes within the boundary 
regions between the domains. At the boundary, the mem- 
branes are bent to connect a TCR/MHCp domain with 
a membrane separation of around 15 nm with an LFA- 
l/ICAM-1 domain with a separation of around 40 nm. 
Similar line tensions, or energies per length, arise at the 
rim of the contact zone, both for TCR/MHCp and LFA- 
l/ICAM-1 domains adjacent to the rim, which we de- 
note here by Xtm and Xli- Of special interest here is 
the difference A r = Xtm — Xli, the energy per length 
for replacing an LFA-l/ICAM-1 boundary at the rim of 
the contact zone by a TCR/MHCp boundary. While A 
is always positive, reflecting the phase separation, A r in 
principle can have both positive and negative sign. 3 
Let us assume that the final TCR/MHCp domain has 



a smaller area than the LFA-l/ICAM-1 domain, which 
seems to be the case for T cells. If A r is larger than A, 
boundaries of TCR/MHCp domains inside the contact 
zone with the LFA-l/ICAM-1 domain are energetically 
more favorable than boundaries at the rim. Hence, the 
final TCR/MHCp domain should be circular and located 
anywhere inside the contact zone, to minimize the overall 
line tension, see Fig. El I n contrast, if A is larger than A r , 
TCR/MHCp domain boundaries at the rim of the contact 
zone are more favorable than interior boundaries, and 
the final TCR/MHCp domain should be in contact with 
the rim. In our simulations, A is clearly larger than A r , 
although A r is positive since our boundary conditions, a 
rim separation of 100 nm, favor LFA-l/ICAM-1 domains 
at the egde of the contact zone. In the case of cells, it 
is reasonable to assume that A is much larger than A r , 
since separation differences between 15 and 40 nm at the 
edge of the contact zone should not cause large energetic 
differences in the cell elasticity. 



3 These line tensions determine the contact angle 9 of the green 
lense in the right pattern of Fig.|H]via cos(0) = [Ali — Xtm)/^- 
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IV. ADHESION DYNAMICS WITH ACTIVE 
TRANSPORT OF TCRS 



In T cells, active processes transport rec eptors into 
the contact zone Wiilfi ng and Davis 19981 and gly- 
coproteins out o f this region [ATIenspach et al. 20011 
iDelon et al. 200lj . The framework enabling these trans- 
port processes is the actin cytoskeleton which polarizes 
during adhesion around the center of the contact zone 
[Alberts et al. 1994 iDustin et al. 1998^ . For TCRs, the 
transport is mediated by myosin, a molecular motor pro- 
tein binding to the actin filaments. Here, we model the 
transport of TCRs as directed diffusion. For simplicity, 
we assume that each TCR molecule experiences a con- 
stant force which is directed towards the center of the 
contact zone midpoint. This force corresponds to an ad- 
ditional term F ■ r in the configurational energy of each 
TCR where F is the magnitude of the force and r the 
distance of the receptor from the center of the contact 
zone. Under the influence of this force, diffusive steps 
bringing TCRs closer to focal point of the cytoskeleton 
in the center of the contact are, in general, more likely 
than diffusive steps in the opposite direction. 

Figure compares the pattern evolution at zero force 
with patterns at the forces F = 0.0lk B T/a ~ 6 • 10~ 16 N 
and F = 0Ak B T/a ~ 6 • 10~ 15 N. The concentrations 
and binding energies are the same as in Fig. [21 For these 
values, the force F = O.OlksT/a is close to the force 
threshold leading to a target-shaped final synapse with 
central TCR/MHCp cluster, see Fig. Besides lead- 
ing to a central TCR/MHCp cluster, the active forces 
speed up the pattern evolution. At the weaker force 
F = O.OlfcsT/a, the final equilibrium state is reached af- 
ter approximately 30 minutes, while the 10-fold stronger 
force F = O.lksT/a leads to equilibrium within few 
minutes. The absolute times are based on the estimate 
that one Monte Carlo step roughly corresponds to 1 ms, 
see above. However, TCR/MHCp rings at intermediate 
times form in all three cases shown in Fig. Only sig- 
nificantly stronger forces prohibit the formation of inter- 
mediate rings by pulling the TCRs more directly to the 
center of the contact zone. 

To quantify the impact of the active forces on the pat- 
tern evolution, we consider again the peripheral ring oc- 
cupation Y of the TCR/MHCp complexes, see Fig. |H1 
Values Y > 80 % at intermediate times indicate a periph- 
eral TCR/MHCp ring as in the inverted T cell synapse, 
and values Y ~ % at later times correspond to target- 
shaped patterns with a central TCR/MHCp domain as in 
the 'mature' T cell synapse. We consider the same molec- 
ular concentrations and binding energies as in the Figs. E| 
and For these values, active forces F < 0. 003 fe^T/a 
are too weak to cause a final central TCR/MHCp clus- 
ter, while forces F > O.lksT/a guide the TCRs very 
quickly to the contact zone center, which prevents the 
formation of a peripheral TCR/MHCp ring at interme- 
diate times. We obtain a T-cell like pattern evolution 
for forces 0.01k B T/a < F < 0.1k B T/a with an inter- 
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FIG. 8: TCR/MHCp ring occupation Y at various forces 
F — f ■ ksT/a as a function of time t for the same molec- 
ular concentrations as in Figs. [3] and [7| A ring occupation 
Y — 100 % corresponds to fully closed peripheral ring of 
TCR/MHCp complexes, smaller percentages of Y correspond 
to partial occupations of the peripheral ring with distances 
r > 35 a from the center of the contact zone. T-cell like 
pattern evolution with intermediate values of Y > 80 % 
(inverted synapse) and final values of Y ~ % ('mature' 
synapse) are obtained for forces between F = O.OlfcsT/a and 
F — 0.1fcsT/a, see also Fig.Q Larger forces prevent the for- 
mation of an intermediate peripheral TCR/MHCp ring, while 
smaller forces do not lead to a final central TCR/MHCp clus- 
ter. The data points represent averages over 24 independent 
Monte Carlo runs for each force. 



mediate inverted synapse and a final 'mature' synapse 
exhibiting a central TCR/MHCp cluster. Experimen- 
tally, the mature synapse of T cells has been observed 
to form on timescales between 5 and 30 minutes, which 
agrees with the equilibration times we obtain for these 
forces. In the absence of active forces (F = 0), the 
intermediate peripheral TCR/MHCp ring seems to be 
metastable and appears in our simulations for times up 
to an hour. This metastability might explain the inverted 
NK cell synapse which consists of a peripheral ring of 
short receptor/ligand complexes, and a central domain 
containing the longer integrins. The inverted synapse 
of NK cells seems to formed by self-assembly, since it is 
not affected by ATP depletion or cyt oskeletal inhibitors 
[Davis et al. 19991 iFassett et al. 2001) . 



V. CONCLUSIONS 

In this article, we have considered possible mechanisms 
underlying the pattern formation during T cell adhesion. 
We propose a novel mechanism for the formation of in- 
termediate patterns, which is based on the nucleation of 
TCR/MHCp microdomains throughout the contact zone 
and the diffusion of free receptors and ligands into this 
zone. This self-assembly mechanism does not require 
active processes and leads to the inverted synapse pat- 
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tern of T cells if the TCR/MHCp concentration is large 
enough. The mechanism leads to multifocal intermedi- 
ates as observed during thymozyte adhesion for smaller 
TCR/MHCp concentrations. 

According to our model, the final T cell pattern with 
a central TCR/MHPp domain is caused by active trans- 
port of TCRs towards the center of the contact zone, 
where the polarized cytoskeleton has a focal point. This 
seems to be in agreement with experiments which show 
that the central TCR/MHCp movement is inhibited by 
blockin g the active cytoskele tal transport with Cytocha- 
lasin D [Grakoui et al. 199Sj . 

In our model, the coalescense of domains eventually 
leads to a single TCR/MHCp domain. Without active 
transport, we find that this domain is located at the rim 
of the contact zone. We obtain the final bull-eye pat- 
tern of the mature synapse with a central TCR/MHCp 
domain only in the presence of active TCR transport. 
For simplicity, we have modeled the active transport by 
a constant average force F which acts on the TCRs and 
is directed towards to the center of the contact zone. We 
obtain the characteristic pattern inversion of T cells for 
average forces in the range 6 ■ 10~ 16 N < F < 6 ■ 10~ 15 N. 
Smaller average forces do not lead to the mature T cell 
synapse with central TCR/MHCp cluster, while larger 
forces disrupt the peripheral TCR/MHCp ring of the in- 
termediate, inverted synapse. It is important to note that 
these average forces acting on a single TCR are signifi- 
cantly smaller than the maximum forces around 1 pN= 
1CP 12 N which can be exerted by a single molecular mo- 
tor. The transport of a TCR molecule over larger dis- 
tances presumably involves several cytoskeletal binding 
and unbinding events. Since the TCRs are membrane- 
anchored, the transport along cytoskeletal fibers close to 
the membrane may also be affected by additional friction 
within the membrane. 

Active cytoskeletal processes could also stabilize the 
multifocal patterns of thymozytes which have been ob- 
served by Hailman et al. The cytoskeleton of thymozytes 
presumably remains in a mobile , nonpolarized stat e 
which still allows cell migration [Hailman et al. 2002j| . 
The few TCR/MHCp clusters of thymozytes may be 
coupled to the cytoskeleton, thus following its move- 
ments. An alternative explanation for the multifocal 
patterns has been given in terms of critic al fluctuations 
[S.-J. Lee et al. 20031 iRavchaudhuri et al| . Fluctuations 
close to a critical point can lead to the appearance and 
disappearance of small domains, since the line tension of 
the domain boundaries then is low. However, the few 
small TCR/MHCp clusters observed bei Hailman et al. 
are rather circular, which indicates a relatively large line 
tension. In addition, the multifocal patterns were ob- 
served over a 100-f old range of antigen concentrations 
|Hailman et al. 2002| . whereas critical fluctuations can 
only be observed in a rather narrow concentration range 
close to the critical point. 

Natural killer (NK) cells form an inverted synapse, 
consisting of a peripheral ring of short HLA-C/KIR 



complexes and a central domain with the longer LFA- 
1/ICAM-l complexes. The formation of the NK cell 
synapse seems not to depend on active cytoskeletal pro- 
cesses, since ATP depletion or disr uption of the cy- 
toskeleton has no e ffect on the pattern [Davis et al. 19991 
iFassett et al. 20 01]. A possible explanation for the NK 
cell synapse is the metastability of the inverted pattern in 
the absence of active cytoskeletal processes. Without ac- 
tive transport, the inverted intermediate synapse persists 
up to an hour in our model. As mentioned above, exper- 
iments show that the central TCR/MHCp movement of 
T cells requires cytoskeletal transport as well. 

We have applied our model here to T cell adhesion, 
using the specific lengths of the TCR/MHCp and ICAM- 
1/LFA-l complexes in the interaction potentials J3J) and 
However, the model is rather general and also applies 
to other cell adhesion events. We have previously consid- 
ered a s impler membrane s ystem with 'stickers' and 're- 
pellers' [Weikl et al. 2002aj| . The phase separation into 
sticker- and repeller-rich domains is driven by the length 
difference between the two molecule types. In the cell 
adhesion geometry, we obtained intermediate patterns 
which are similar to those presented here. A difference 
to T cell membranes is that the repeller-rich domains are 
unbound. Large-scale membrane fluctuations in these 
domains then drive the final sticker clusters towards the 
center of the contact zone, at least for the 'free' bound- 
ary conditions with unc onstrained membra ne separation 
at the contact zone rim [Weikl et al. 2002a| . In contrast, 
the coexisting TCR/MHCp and ICAM-l/LFA-1 domain 
types of T cells are both bound, and large-scale mem- 
brane fluctuations are suppressed. 

We have characterized the receptor/ligands by an ef- 
fective binding energy which is directly related to an ef- 
fective 2D equilibrium constant Kid, see Appendix. The 
equilibrium constant Kid is the ratio of the kinetic on- 
and off-rates k on and fc„/ / for the receptor-ligand bind- 
ing [Bell. 1 978], Characterising the binding kinetics by 
the single parameter K%d rather than the two parame- 
ters k on and fc // is justified at least if the on-reaction 
of a receptor-ligand pair in apposing membrane patches 
of our discrete model is faster than the timescale 1 ms 
for the diffusive Monte Carlo steps. In general, such an 
approach may also be justified by a local equilibration 
within domains. 



APPENDIX A: BINDING ENERGIES AND 2D 
EQUILIBRIUM CONSTANTS 

In this appendix, we consider the relation between the 
2d equilibrium constants and the binding energies of re- 
ceptors and ligands with square- well potentials J2J| or 10}. 
The 2d equilibrium constants are defined by 



xrxl 

where xr is the area concentration of receptors R in one 
of the membranes, xl is the area concentration of the 
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ligands L in the apposing membrane, and x R l is the con- 
centration of the complexes. In general, K 2 d depends on 
the state of the membrane, not only on the interaction 
potential of receptors and ligands. For example, if the 
membrane separation is to large to allow complex forma- 
tion, the equilibrium constant (|A1|> is zero. In a bound 
state, K<id will depend on the separation and roughness 
of the membrane, which in turn are affected by the con- 
centrations of the receptors, ligands, and steric repellers 
such as glycoproteins. Here, we only consider the 'ideal' 
2d equilibrium constant for a membrane segment which is 
entirely within the binding range of the receptor/ligand 
square- well interaction. 

For convenience, we choose in the following the grand- 
canonical ensemble with chemical potentials \i R and \il 
for receptors and ligands. However, the ideal 2d equilib- 
rium constant derived below will be independent of [i R 
and hl and thus applies also to the canonical ensemble 
with fixed overall receptor and ligand concentrations. 

Let us consider two apposing membrane patches within 
binding range of the receptors R and ligands L. A 'state' 
of these apposing patches is characterized by the numbers 
niL and m R of ligands and receptors present in the two 
membranes. The configurational energy is h{mL,m R ) = 
Url min(m L , m R ) - [ihm L - fi R m R where U RL is the 
binding energy of RL complexes. To simplify the nota- 
tion below, the parameters fi R , fiL, and Url are taken 
to be in units of the thermal energy k R T. In a given 
state, there can be k bound RL complexes, plus cither i 
uncomplexed ligands or j uncomplexed receptors (or no 
additional uncomplexed molecules). The partition func- 
tion then has the form 



oo oo 



= ^ ^2 cx P(- h ( m L,m R )) 

m,L—0 mji—0 

CO / OO OO \ 

_ e HpL+^R-u RL ) 1 1 -|- e i/J,L + e 3 J 

i=l j=l ' 

^ - A (A2) 



fc=0 



X — gUL+^R — URL \ 1 — e^-L \ _ gMfi 



and the concentration of bound RL complexes is 



XRL = -L >T fcek^+fK-tf**) ( 1 + £ e 



k=l 



i=l 

oHl+Hr — Url 



(A3) 



In deriving these expressions we made use of 



k=0 k=0 



kgk = (A4) 



for s < 1. The concentration of uncomplexed receptors 



is given by 



xr 



oo oo 



k=0 

-(— 

a 2 \ 1 - e» R 



1 



\ — e ^R+fJ.L 



(A5) 



and accordingly, the concentration of uncomplexed lig- 
ands is 



1 

XL = o 



1 



1 



(A6) 



Hence, the ideal 2-dimensional equilibrium constant of 
receptors and ligands in apposing membrane patches 
within binding range is given by 



K- 



2d 



X R L 



(1 - e^+t 1 *) 



a 2 e- u *i 



(A7) 



The last expression holds for e flL+fJ ' R <C 1 which is true 
for the receptor and ligand concentrations studied in this 
article, see below. 

For comparison, let us also consider two apposing 
membrane patches with a separation outside of the bind- 
ing range of receptors and ligands. The partition function 
for the patches then is 



1 



i=0 



3=0 



(1 - e^)(l - e^«) 



(A8) 



and the concentrations of receptors and ligands are given 
by 



x R 



Eoo 



XL 



a 2 (1 - e^) 
(A9) 
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